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Digital pelvic radiographs are used to identify the locations of implanted iodine- 125 seeds and their 
numbers after insertion. However, it is difficult and laborious to visually identify and count all implanted 
seeds on the pelvic radiographs within a short time. Therefore, our purpose in this research was to develop 
an automated method for estimation of the number of implanted seeds based on two-view analysis of pelvic 
radiographs. First, the images of the seed candidates on the pelvic image were enhanced using a difference 
of Gaussian filter, and were identified by binarizing the enhanced image with a threshold value determined 
by multiple-gray level thresholding. Second, a simple rule-base method using ten image features was 
applied for false positive removal. Third, the candidates for the likely number of a multiply overlapping 
seed region, which may include one or more seeds, were estimated by a seed area histogram analysis and 
calculation of the probability of the likely number of overlapping seeds. As a result, the proposed method 
detected 99.9% of implanted seeds with 0.71 false positives per image on average in a test for training 
cases, and 99.2% with 0.32 false positives in a validation test. Moreover, the number of implanted seeds 
was estimated correctly at an overall recognition rate of 100% in the validation test using the proposed 
method. Therefore, the verification time for the number of implanted seeds could be reduced by the provi- 
sion of several candidates for the likely number of seeds. 

Keywords: prostate brachytherapy; estimation of number of iodine-125 seeds; pelvic radiographs; image 
processing in radiotherapy; two-view analysis 



INTRODUCTION 

The incidence of prostate cancer is rising rapidly not only 
in western countries but also in Asian countries, including 
Japan [1, 2]. About 32 000 deaths from prostate cancer 
were anticipated in the USA in 2010 [3], and about 10 000 
deaths in Japan in 2007. A number of patients with early 
prostate cancer are treated by means of brachytherapy with 
permanent seed implants including iodine-125 seeds at low 



dose rates [4]. The advantage of this implant procedure is 
that there is sufficient dose coverage for localized prostate 
cancer, and reduction of harmful effects on healthy tissue. 
Furthermore, brachytherapy with iodine-125 seeds can pre- 
serve sexual function compared with external radiation and 
surgical prostatectomy [5-7]. 

The brachytherapy procedure with iodine-125 seeds con- 
sists of three steps: planning, implant of seeds and post- 
operative evaluation. The postoperative evaluation after 
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seed implantation is critical for validating the correct place- 
ment and the number of the implanted seeds according to 
preoperative or intraoperative planning, because seeds can 
migrate due to edema and swelling in the days immediately 
following their insertion [8, 9]. Seed migration could in- 
crease the probability of recurrence because of inadequate 
dose distribution in the prostate [10-13], and this could 
cause complications such as pulmonary and coronary artery 
embolization [14-18]. In this study, we did not deal with 
verification of correct placement of iodine-125 seeds; 
rather, this study focuses on estimation of the number of 
implanted seeds. The postoperative evaluation validating 
the number of the implanted seeds should be performed 
within a short time for the following two reasons: (i) a dis- 
agreement in the number of seeds between radiographs 
would be treated as an incidence of radiation source loss 
according to the Japanese Laws Concerning the Prevention 
of Radiation Hazards due to Radioisotopes and and; (ii)any 
adverse patient reactions to seed implantation, such as 
blood pressure reductions, body temperature decrease and 
hypopnea can be dealt with straight away in a hospital 
environment. 

If the number of seeds on post-insertion pelvic radio- 
graphs of a patient does not agree with the true number of 
seeds, the patient cannot be moved out of the operating 
room. Therefore, the location and number of implanted 
seeds should be correctly and instantly identified just after 
seed insertion prior to moving the patient out of the operat- 
ing room. 

For this purpose, a number of automated identification 
methods for locating radioactive iodine-125 seeds and 
ascertaining the number of seeds for prostate brachytherapy 
have been developed based on three-dimensional (3D) 
reconstructed images by use of fluoroscopy images or digi- 
tized radiographs [19-21], computed tomography (CT) 
[22-24], magnetic resonance imaging (MRI) [25, 26] and 
tomosynthesis [27, 28]. However, it is difficult to employ 
imaging modalities such as CT or MRI in the operating 
room. Furthermore, since the 3D images are constructed 
using more than two views, obtained at different times, 
prostate displacement of the seeds due to motion and respir- 
ation can mean that they cannot be correctly detected [29, 
30]. For that reason, two-view pelvic digital radiography is 
commonly used in clinical situations in Japan. However, it 
is difficult and laborious to visually and manually identify a 
large total number of implanted seeds (e.g. -70) in digital 
radiographs within a short time span of 1 or 2 min, because 
the seeds are very small (around 1 mm diameter x 5 mm 
length), and two or more seeds often overlap in a projection 
image. Therefore, further studies on methods for automatic- 
ally and immediately verifying the number of implanted 
seeds based on a simple radiograph are needed for assisting 
radiation oncologists or medical physicists in the post- 
operative evaluation procedure. 



Our purpose in this research was to develop a fully auto- 
mated method for estimation of the number of implanted 
iodine-125 seeds instantly prior to moving a patient out of 
the operating room, based on two-view analysis of pelvic 
radiographs, obtained just after seed insertion. 

MATERIALS AND METHODS 
Clinical cases 

This study was approved by the Institutional Review Board 
(IRB) of our university hospital. We started to treat early- 
stage prostate cancer patients using iodine-125 brachyther- 
apy from 2007. Twenty-eight patients (ages: 53-74 years, 
mean: 65 years) were selected from all patients from 2007 
to 2009 in our university hospital, and this duration was 
approved by the IRB. All patients were treated with 
iodine-125 seeds (0.97 mm diameter x 4.55 mm length, 
OncoSeed GE Healthcare, Inc., Arlington Heights, USA). 
Fifty-six pelvic digital images of the 28 patients were 
acquired after seed implantation by a digital radiography 
(DR) system (Konica Minolta Health Care Co., Ltd, Tokyo, 
Japan) in two projections by use of an X-ray mobile C-arm 
(Toshiba Medical Systems Co., Ltd, Tokyo, Japan). The in- 
cidental X-ray angles were decided within an arbitrary ± 
15° so there would be as little seed overlap as possible with 
the use of fluoroscopy and so as not to touch the patient or 
surgical table. The pelvic images were taken with an anti- 
scatter grid (grid ratio 5:1, grid density 34 lines/cm) at 
X-ray tube voltages of 80-90 kV and 10-20 mAsecond. 
The matrix size, pixel size and the number of gray levels in 
the pelvic DR images were 1430x1722, 0.175 mm, and 
4096, respectively. Forty-nine pelvic radiographs were 
selected after excluding seven images, which had a lot of 
noise, blurred edges due to patient movement, and 
implanted seeds outside the irradiation field for this study. 
A total of3582 of the implanted iodine-125 seeds were 
visually confirmed in the 49 pelvic radiographs with con- 
sensus by two radiation oncologists and a medical physi- 
cist. The mean number of seeds was 73 per patient. For a 
validation test of our proposed method, we divided the 28 
cases including 49 images into two groups, i.e. training and 
test images; the numbers of images were 21 and 28, 
respectively. 

Automated identification of implanted seeds 

Figure 1 shows a schematic diagram of the automated iden- 
tification of implanted seeds. The seed candidates were 
identified using the following steps. First, a field of view 
with a circular shape was segmented by analysis of edge 
pixel-value profiles in the horizontal and vertical directions 
within the edge-enhanced image by a Sobel filter (mask 
size: 7x7 pixels). The Sobel filters in the horizontal and 
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Fig. 1. Proposed scheme for automated identification of implanted seeds. 
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The large mask size oil xl was employed for detection 
of large structure edges such as the irradiation field. Only 
magnitudes of the gradient vectors obtained by the Sobel 
filters were used in enhancement of the image edges. 
Second, anti-scatter grid lines on the image were suppressed 
by a low pass filter (below 0.56 mm -1 ). Third, the images 
of the seeds were enhanced using a difference of Gaussian 
(DOG) filter, i.e. subtraction of a smoothed image with a 
Gaussian filter (mask size: 25 x 25 pixels, standard devi- 
ation: 17 pixels) from the suppressed grid-line image. The 
large Gaussian filter with a mask size of 25 x 25 pixels was 
used to suppress the implanted seeds as much as possible. 
Finally, numerous seed candidate regions were obtained by 



binarizing the enhanced seed image with a threshold value 
determined by use of a multiple-gray-level thresholding 
technique. The threshold value was determined so that the 
number of seed candidates would be smaller than 160. On 
the other hand, the edge-enhanced image with a Sobel filter 
(mask size: 7x7 pixels) was binarized by one of the three 
threshold values, which was automatically determined by 
one of the ratios of 0.06, 0.015 and 0.005 peak level in the 
pixel histogram based on the threshold value of the 
enhanced image of the seed. Each peak level was deter- 
mined empirically in this study. The grid lines were sup- 
pressed by the low pass filter below 0.56 mm -1 . Figure 2 
shows illustrations for segmentation of seed candidate 
regions based on the DOG filter and a multiple-grey level 
thresholding technique. By applying the DOG filter, i.e. 
subtraction of the Gaussian-smoothed image (Fig. 2b) from 
the grid line-suppressed image (Fig. 2a), background struc- 
tures such as the pelvic bone, colon gas and the bladder 
were attenuated effectively with contrast medium, as shown 
in Fig. 2c, whereas the seed contrast became higher than 
that in the original image. Finally, as shown in Fig. 2d, seed 
candidate regions were segmented well, including some 
obvious false positive (FP) regions, by use of the 
multiple-gray-level thresholding technique. The binarized 
edge-enhanced images were used only for determination of 
region coincidence, which is one of the image features for 
FP removal. 

False positive removal by a rule-based method 

A simple rule-based method was applied for removal of 
seed FPs with use of nine image features, i.e. minimum, 
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Fig. 2. Illustrations of segmentation of seed candidate regions: (a) a grid-line-suppressed image, 
(b) an image smoothed by a Gaussian filter, (c) an enhanced image of a seed, (d) seed candidate regions 
obtained by a multiple-gray-level thresholding technique. 



maximum and mean pixel values and their SD, contrast in 
the enhanced seed image, area (the number of pixels), circu- 
larity, slenderness and region coincidence. Contrast was 
defined as the maximum pixel value minus the minimum 
pixel value in the seed candidate region. Circularity was 
defined as the ratio of the overlapping area of a seed candi- 
date within the equivalent-area circle to the area of the seed 
candidate. Region coincidence was defined as the ratio of 
the logical AND region to the logical OR region between 
a binarized seed enhancement image and a binarized 
edge-enhancement image. Each rule for each image feature 
used in the rule-based scheme was based on a simple thresh- 
old value, which was determined from all true seed 
regions. The rules for the minimum and maximum pixel 
value were greater than each mean value plus 2 SD and 3 
SD, respectively, and the rules for the mean pixel value, SD 
and contrast were less than each mean value minus 3 SD. 
Seed candidates with a number of pixels >50 and <1400 



were removed as noise or large structures of pelvic bone or 
X-ray marker, respectively. The rule of the circularity for FP 
removal was below 10% or above 90%. Slenderness was 
defined as the number of times morphological erosion of a 
binary image was performed until the seed candidate disap- 
peared, and the rule for FP removal was a value less than 
the mean minus 2 SD. The region coincidence rale for FP 
removal was <40%. 

Outlier FP removal by accumulated profile analysis 

The remaining outlier FPs on the bone edge, which exist 
outside a 'prostate' region, were removed by use of accu- 
mulated profiles based on the number of pixels in a binary 
image with segmented seed candidate regions in the hori- 
zontal and vertical directions. The 'prostate' region is 
where a prostate is supposed to be. Figure 3 shows accumu- 
lated profiles on the binary image in two directions. The 
bin width for producing the accumulated profile was set as 
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Fig. 3. Accumulated profiles on a binary image with segmented 
seed candidate regions in horizontal and vertical directions. 

70 pixels (12.25 mm) because the seeds were inserted at 
10-mm intervals with use of a perineal template, and the 
magnification ratio of the seeds on an image was almost 
1.2x in this radiographic geometry. The maximum value of 
the accumulated profile of each direction was detected, and 
then a circumscribed quadrangle around the 'prostate' was 
determined by four positions corresponding to four feet of 
two profile peaks <5% of the maximum value in the hori- 
zontal and vertical directions, as shown in Fig. 3. Finally, 
several outlier FPs outside the 'prostate' region were 
removed by use of the circumscribed quadrangle. 

Performance evaluation of automated identification 
of implanted seeds 

The center coordinates of all seeds on a pelvic image were 
determined by two radiation oncologists and a medical 
physicist, and considered as gold standards. A seed candi- 
date was detected as true positive (TP) if the seed candidate 
region included the 'center' coordinate, but otherwise was 
considered as a false negative (FN). 



Estimation of likely numbers of 
MOS regions 

1 

Estimation of probability of each 
number of MOS regions 

r 

Several candidates of the most likely number 
of MOS regions 

Fig. 4. Proposed scheme for automated estimation of the number 
of seed regions including isolated seed regions and multiply 
overlapping seed (MOS) regions. 

of the variation of pixel values on the periphery 
and overlapping areas in the MOS region. 

(iii) The likely numbers of MOS regions including 
two or more seeds were calculated. 

(iv) The probability of each number of overlapping 
seeds was estimated based on the frequency of 
the likely numbers of the seeds. 

(v) The most likely candidates for the number of 
overlapping seeds were selected based on their 
higher probabilities in one-view analysis of pelvic 
radiograph. 

(vi) The most likely candidates for the number of 
overlapping seeds were selected based on their 
higher probabilities in two-view analysis of 
pelvic radiographs. 



Automated estimation of the number of implanted 

seeds 

Figure 4 shows a proposed scheme for automated estima- 
tion of the number of implanted seeds. The candidates for 
the likely number of each region, which may include one 
or more seeds, were determined based on a seed candidate 
area histogram using the following six steps. Note that only 
true positive seeds were dealt with in this step. 

(i) The number of isolated seeds was estimated from 
a major peak in a seed candidate area histogram. 

(ii) The number of pixels in each multiply overlapping 
seed (MOS) region was corrected by consideration 



Estimation of the number of isolated seeds based 
on seed candidate area histograms 

Figure 5 shows an example of the seed candidate area 
histogram for all cases used in this study. The seed candi- 
date area histogram was used for classifying all seed candi- 
dates into isolated seeds and overlapping seeds, and for 
calculating the mean and SD of isolated seeds area. We 
assumed that the seed regions smaller than a threshold area 
at the largest peak of the seed candidate area histogram 
were considered as isolated seeds. That is because original- 
ly all seeds were individually implanted within a prostate, 
and thus the largest peak is considered to include a majority 
of isolated seeds. 
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Correction of the number of pixels in a region 
withMOS 

Basically, the number of seeds in an MOS region is esti- 
mated by dividing a predicted number of pixels within the 
MOS region by the mean area of the isolated seed regions. 
However, there are two requirements for more accurate esti- 
mation of the number of seeds in the MOS region if the 
number is estimated in one of the projection images as 
follows: Requirement (i) each MOS region should be reseg- 
mented because the pixel values on the periphery of the 
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170 220 270 320 
Area (number of pixels) 

Fig. 5. Seed candidate area histogram for all cases used in this 
study. The seed regions smaller than a threshold area at the largest 
peak of this histogram were considered as isolated seeds. 



MOS region vary with the seed location, tube voltage or 
abdominal thickness, and thus the contrast of each MOS 
region to the background in a projection image could 
change with such factors; and (ii) Requirement (ii) we need 
to count the number of pixels in an overlapping area twice 
by identifying the overlapping area in one projection of the 
pelvic image with implanted seeds. For example, if the 
number of pixels in a region with two overlapping seeds 
was estimated in a projection image, the number of pixels 
in the overlapping seed region would be smaller than twice 
the mean area of the isolated seed regions due to an over- 
lapping area between the two seeds. 

Therefore, the number of pixels in the MOS regions 
obtained by the multiple-gray level thresholding explained 
in the section on 'Identification of initial seed candidate 
regions' was corrected using the following method. 
Figure 6 shows the correction of the number of pixels in an 
MOS region based on the following two steps. In the first 
step for Requirement (i), a threshold value was determined 
for each MOS region by changing the mean pixel value of 
the isolated seeds plus |3 x SD, where (3 was changed from 
0.1 to 5.5 in increments of 0.1, until the average incremen- 
tal number of pixels for all MOS regions was equal to the 
average number of pixels on the peripheral lengths of all 
isolated seeds. Then, all MOS regions were resegmented by 
the threshold value as shown in Figure 6c. In the second 
step for Requirement (ii), the overlapping areas were identi- 
fied based on a thinned image for each MOS region. Let us 
consider the case where a seed overlaps with another seed. 
In this instance, the pixel values on the centerline, which 
correspond to an overlapping area between the two seeds, 
are slightly lower than those on the centerline where they 
are not overlapping in an enhanced seed image. Therefore, 




(e) 



(f) 



□ initially segmented region 

H overlapped area 

■ grown MOS region in the first step 



Fig. 6. Correction of the number of pixels in a multiply overlapped seed (MOS) region: (a) an enhanced image of a 
seed with an MOS region including five seeds, (b) an initially segmented MOS region derived from multiple-gray 
level thresholding, (c) a grown MOS region in the first step, (d) a thinned image derived from the initial MOS region 
(b), (e) a segmented overlapped area in the second step and (f) an image including the results of the two steps. 



748 



H. Kawata et al. 



a thinned image was obtained for determination of center- 
lines from the initially segmented MOS region by 
multiple-gray-level thresholding. Next, pixels, which are 
lower than the mean pixel value minus the SD on the cen- 
terlines, are considered as pixels in overlapping areas. 
Finally, the number of such pixels is counted twice. 

Figure 6a shows an enhanced seed image with an MOS 
region including five seeds. Some areas, which may be 
overlapping areas, are darker than other areas. Figure 6d 
shows a thinned image derived from the initial MOS 
region, shown in Fig. 6b. Figures 6e and f show a segmen- 
ted overlapping area in the second step and an image in- 
cluding all results in the two steps. 



Estimation of the likely number of overlapping 
seeds based on one-view analysis of pelvic 
radiograph 

The probability P(x) of likely numbers x of overlapping 
seeds including two or more was defined by the frequency 
of the number of seeds determined from the candidate area 
divided by the mean area of isolated seeds with changing 
SDs. The probability was calculated for each MOS region 
using the following algorithm: 

Step 1: Setting a and Nc as -4.0 and 0, respectively, 
and calculation of the average area of the isolated 
seeds. 

Step 2: Increment of a and Nc by 0.1 and 1, 
respectively. 

Step 3: Calculation of the likely number x of overlap- 
ping seeds using the following equation: 

A c 

x = ^ 

A + aa 

where [ ] is the floor function, A c is the candidate 
area and Aj is the isolated candidate area. 

Step 4: Count of the frequency Nx of each likely 
number of overlapping seeds. 

Step 5: If rx>4.0, the probabilities of all candidates 
of the likely number x of overlapping seeds are cal- 
culated by Nx/Nc. 

If not, go to Step 2. 

Several candidates for the most likely number of the 
overlapping seeds were selected based on their higher prob- 
abilities. Figure 7 shows an example histogram of the likely 
numbers of an MOS region including five overlapping 
seeds. To evaluate this estimation method of the likely 
number of overlapping seeds, the recognition rate was cal- 
culated by dividing the estimated number of seeds by the 
true number of implanted seeds, which was determined by 
two radiation oncologists and a medical physicist. 
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Likely number of seeds 

Fig. 7. Example of a histogram of the likely numbers of an 
MOS region including five overlapping seeds. 

Estimation of the likely number of overlapping seeds 
based on two-view analysis of pelvic radiographs 

The final numbers of overlapping seeds were estimated 
based on two-view analysis of pelvic radiographs. Because 
the regions of all overlapping seeds on two pelvic radio- 
graphs did not agree with each other, the location and 
number of regions of overlapping seeds are inconsistent 
with each other. Therefore, we adopted two-view analysis 
of pelvic radiographs based on two results of one-view ana- 
lysis of pelvic radiographs in two directions. The most 
likely five candidates for the number of overlapping seeds 
were selected from all candidates of the most likely number 
of the overlapping seeds, which were determined based on 
the one-view analysis of the two pelvic radiographs of a 
patient. With this method, six patients were employed as a 
training data set and a test data set, respectively. 

The candidate for total number of all implanted seeds 
was detected as a TP if the candidate, which was selected 
by its higher probability, included the true number of 
implanted seeds in either result of the two-view analysis of 
pelvic radiographs, but otherwise was an FN. The overall 
recognition rates of the number of implanted seeds were 
calculated by dividing the number of TPs by the number of 
patients. 

RESULTS 

The evaluation was based on two-step evaluations, i.e. (i) 
the sensitivity of seed locations and (ii) the recognition rate 
of the isolated seeds and the likely number of overlapping 
seeds for TP seeds. All procedures took 224.1 ±30.3 s to 
carry out using a central processing unit (CPU) on a 
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personal computer with a 2.53-GHz-CPU (Intel Core 2 
Duo) 4-GB random access memory (RAM). The proposed 
method detected 99.9% (1508/1510) of implanted seeds 
with 0.71 FP per image on average in a test for training 
cases. Table 1 shows the results of automated identification 
of implanted seeds in the training cases and the test cases. 
Twenty-four false positives were detected on the pubic 
bones and urethra catheters, and eighteen seeds in 
very-low-contrast regions were missed. All isolated seed 
candidates (544) were correctly identified as one seed. 



Figure 8 shows the relationship between the number of can- 
didates for the most likely number of seed regions and the 
recognition rate of the automated estimation in the training 
cases and the test cases. The recognition rate of overlapping 
seeds was 72.4% when we used the first candidate of the 
likely number of overlapping seeds in the training cases. If 
the most likely three candidates of the number of overlap- 
ping seeds were selected, the correct number of seeds for 
all regions was included in the three candidates with a rec- 
ognition rate of 97.5% (Fig. 8a). 



Table 1. Results of automated identification of implanted seeds in the training cases and the test cases 
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N umber of candidates selected Number of candidates selected 

in order of likely number of seeds in order of likely number of seeds 

Fig. 8. Relationship between the number of candidates for the most likely number of seed regions and the recognition rate of automated 
estimation: (a) training cases, and (b) test cases. 



In the test cases, the sensitivity of identification of 
implanted seeds was 99.2% (1982/2073) with 0.32 FP per 
test image as shown in Table 1. All isolated seeds (515) 
were identified correctly. The recognition rate for all 
regions was 92.4% using only the first candidate for the 
likely number of overlapping seeds, and 99.1% using three 
candidates for the likely number. 

Figure 9 shows the relationship between the number of 
candidates for the most likely number of all implanted 
seeds and the overall recognition rate based on two-view 
analysis of pelvic radiographs. The overall recognition rate 
was up to 83.3% at the fourth or fifth candidate for number 
of all implanted seeds in the training cases, and was 
100.0% for the third candidate in the test cases. 

DISCUSSION 

Tubic et al. [19] developed an automated method for the 
detection of the position and the orientation of implanted 
seeds on fluoroscopic images or scanned radiographs. They 
reported that the sensitivity for detection of implanted 
seeds was 92.0% using their method. Su et al. [20] devel- 
oped a method of determination of 3D seed positions from 
multiple radiographic projections by dealing with overlap 
of seed images. Their method produced 3.0 FPs with a sen- 
sitivity of 98% in a phantom study (implanted 64 seeds), 
based on 3D reconstruction of three projections of fluoro- 
scopic images. Benkhoucha et al. [28] developed an auto- 
matic tomosynthesis-based seed detection technique, which 
yielded seed a detection rate of 96.7% and 2.6% FPs. 
Whitehead et al. [24] proposed a technique for localizing 



-•■ " Training cases (n=6) B Test cases (n=6> 




Number of candidates selected in order of 
likely number of seeds 

Fig. 9. Relationship between the number of candidates for the 
most likely number of seed regions and the overall recognition 
rate with two-projection pelvic images analysis. 



brachytherapy seeds by utilizing voxel intensity, which 
would provide seed detection and localization especially 
for low resolution CT data sets. In four separate simulations 
(X-ray CT images of three separate phantoms containing 
multiple implanted seeds in different configurations, and in 
an in vivo X-ray CT volume data set), all simulated seed 
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objects were detected. On the other hand, our proposed 
method detected 99.9% of implanted seeds with 0.71 FPs 
per image on average in a test for training cases (21 clinical 
images), and 99.2% with 0.32 FPs in a validation test (28 
clinical images). Therefore, the detection sensitivity seems 
to be comparable with the past studies, and the number of 
false positives is improved. 

Eighteen seeds (0.5%), which were overlapping with 
anatomic regions or materials such as pubic bone or a cath- 
eter, were missed by the automated identification of 
implanted seeds in the proposed method. Especially, in a 
radiograph of one particular patient (#12), 12 FNs occurred 
in lower contrast regions where the pubic symphysis and 
the bladder with contrast media were overlapping with 
implanted seeds. In addition, the edges of the pelvic bone, 
as well as parts of the pubic symphysis and pelvic cavity, 
were detected as 0.32 FPs per image because of their simi- 
larity to implanted seeds. Therefore, we should improve the 
automated seed identification method in terms of the detec- 
tion of seeds and removal of FPs in future work. 

The overall recognition rate of number of implanted 
seeds in the test cases was 100% using two- view analysis 
of pelvic radiographs, if three candidates were selected for 
the likely number. Therefore, it is possible to use this tech- 
nique for final verification of the number of implanted 
seeds. Furthermore, the proposed method may reduce the 
radiologists' workloads in a clinical situation, as radiolo- 
gists would only have to count the number of isolated 
seeds, but they should verify the numbers of multiply over- 
lapping seeds. The overall recognition rate could be 
improved if using several projection images (three or more) 
taken at various angles [21], but patient dose and labor for 
radiography would be increased compared with this study. 
Our proposed method, based on two-view analysis of 
pelvic radiographs, gave a good estimate of the number of 
implanted seeds with minimal patient dose and labor for 
radiography. 

As shown in Fig. 8b, all candidates for the most three 
likely numbers estimated by the proposed method were in- 
correct for 0.9% of overlapping seeds in the test cases. 

In this study, three candidate numbers for each MOS 
region with more than one seed were selected based on the 
three highest probabilities. This depends on the segmentation 
accuracies of overlapping seed regions and isolated seed 
regions. Furthermore, some of the overlapping seed regions 
were not segmented well, and the most likely candidate 
numbers tended to be incorrect. Therefore, we should 
improve the segmentation method for the overlapping seed 
regions in future work. 

We recognize that there are two limitations in the pro- 
posed method. First, the users of the proposed method have 
to select one of no more than three candidate numbers for 
the multiply overlapping seeds produced by the proposed 
method. This would not be practical in a clinical situation. 



One possibility to overcome this would be to verify whether 
the 'true' number of seeds, which the radiologist implanted, 
is included in the list of candidate numbers of implanted 
seeds suggested by the proposed method. Second, there is 
an average computation time of about 4 min, which may be 
relatively long for clinical use, but as the time could be 
greater shortened to <30 s if graphical processing units 
(GPUs) were to be employed. In addition, we should 
improve the algorithm of the proposed method with respect 
to efficient computation. 

CONCLUSIONS 

Because a disagreement in the number of seeds between 
radiographs would be treated as an incidence of radiation 
source loss according to the Japanese Laws Concerning the 
Prevention of Radiation Hazards due to Radioisotopes and 
Others, we need an automated method for verifying the 
number of seeds within a short time. Therefore, we have 
developed an automated method for estimation of the 
numbers of implanted iodine- 125 seeds on a pelvic radio- 
graph after permanent prostate implants. The proposed 
method detected 99.2% of implanted seeds with 0.32 false 
positives in a validation test. If the most likely three candi- 
dates for the number of implanted seeds were selected, the 
number of implanted seeds was estimated correctly at an 
overall recognition rate of 100% in the validation test. 
Therefore, the preliminary results showed that the proposed 
method would be useful for estimation of the numbers of 
implanted seeds on pelvic radiographs with a high rate of 
seed detection and recognition of the number of seeds. 
Furthermore, the proposed method would be simple to 
apply in clinical situations Justas it only requires two-view 
pelvic radiographs. 
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